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If  a  steel  rod  is  hung  by  a  thread  and  struck  at  the  lower 
end,  an  audible  metallic  "ring"  is  produced  whose  fundamental  fre¬ 
quency  depends  on  the  velocity  of  sound  in  steel  and  on  the  length 
of  the  rod.  At  the  very  point  of  impact  on  the  rod,  a  compressional 
wave  is  started  which  emenates  outward  and  returns  from  free  surfaces 
as  rarefaction  waves.  Most  of  the  waves  cancel  out  except  those 
whose  wavelength  are  multiples  of  the  rod  dimensions.  The  initial 
energy  imparted  to  the  rod  is  ultimately  dissipated  as  heat  but  only 
after  a  relatively  long  period  of  time.  A  small  amount  of  the  energy 
is  transferred  to  the  air  as  a  sound  (this  is  the  portion  ve 

hear) .  The  waves  we  have  induced  in  the  rod  are  actually  identical 
to  the  very  long  wavelength  phonons  (or  normal  modes)  already  present 
in  the  rod  as  part  of  the  nornal  thermal  spectrum  of  phonons.  When 
the  rod  is  struck,  this  portion  of  the  phonon  spectrum  is  considerably 
enhanced,  but  only  after  many  seconds  does  this  energy  redistribute 
itself  so  that  thermal  equilibrium  is  re-established.  On  the  other 
hand  if  an  X-ray  beam  strikes  the  rod,  it  will  produce  phonons  whose 
frequency  is  M.0i2-1013  hertz,  at  least  109  times  greater  than  the 
audible  frequencies  created  by  mechanical  striking.  Furthermore, 
rather  than  taking  seconds  to  reach  thermal  equilibrium,  the  time 
required  is  only  ■xlO”12  seconds  for  these  high  frequency  phonons. 

The  reason  for  this  difference  in  behavior  between  the  low  frequency 
phonons  and  the  high  frequency  phonons  is  well  understood  in  terms 
of  anharmonicity. 

The  electrons  binding  atoms  together  in  a  solid  can  (in 
a  crude  first  approximation)  be  considered  to  be  replaced  by  springs 
joining  nearest  and  next  nearest-  neighbors.  For  small  atomic  displace- 
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ments  the  forces  created  by  these  springs  are  harmonic  (quadratic 
potential  energy).  Within  this  approximation,  the  atoms  and  springs 
become  a  system  of  coupled  harmonic  oscillators  and  this  gives  rise 
to  the  normal  modes,  (or  phonons)  which  are  excited  according  to 
the  usual  exponential  Boltzman  partition  function  depending  on  the 
temperature  of  the  system  and  the  frequency  of  the  phonon.  However, 
in  reality,  a  system  in  which  the  forces  between  atoms  are  truly 
harmonic  can  never  reach  thermal  equilibrium  since  the  phonons  have 
no  means  for  exchanging  energy  with  each  other.  The  introduction 
of  anharmonicity,  i.e.  higher  energy  cubic  and  quartic  terms  in 
the  potential  energy  of  the  compressed  or  elongated  springs,  provides 
a  mechanism  for  coupling  the  phonons  together.  Such  higher  order 
terms  are  generally  small  but  they  alone  are  responsible  for  such 
physical  properties  as  thermal  expansion,  the  approach  to  thermal 
equilibrium,  and  thermal  conductivity.  The  influence  of  the  cubic 
and  quartic  terms  in  the  potential  energy  depends  on  the  amplitude 
of  the  displacements  being  almcat  negligible  for  small  displacements. 
This  basically  is  why  the  sound  waves  produced  in  the  steel  bar 
ring  for  so  long  a  period  of  time  and  travel  virtually  unattenuated 
at  the  speed  of  sound  over  the  length  of  the  bar.  Associated  with 
these  sound  waves  are  atomic  displacements  less  than  10  of  the 
distance  between,  atoms.  On  the  other  hand  the  amnlitude  of  vibration 
for  thermal  waves  are  at  least  106  orders  of  magnitude  greater  so 
that  anharmonicity  quickly  couples  the  phonons  together.  This  effect¬ 
ively  scatters  these  waves  so  that  the  transfer  of  heat  energy  pro¬ 
ceeds  slowly  via  a  diffusion  process. 

If  our  steel  rod  is  replaced  by  a  TiNi  alloy  called  Nitinol 
and  struck  at  about  100°F  a  similar  metallic  "ring"  is  heard  but 
by  lowering  the  temperature  to  about  80°F  the  rod  can  no  longer 
be  caused  to  "ring".  The  effect  over  such  a  small  temperature  range 
is  quite  dramatic  and  in  fact  even  at  very  high  stress  levels  the 
compressional  waves  are  quickly  damped  and  converted  into  heat. 

The  prospects  of  employing  this  high  damping  capacity  to  absorb 
vibrations  in  mechanical  systems  appear  promising. 

The  reason  for  the  very  effective  damping  in  NiTi  at  room 
temperature  is  related  to  the  crystallographic  (martensitic)  phase 
change  that  occurs.  As  a  precursor  to  this  transition  we  find  a  so- 
called  soft  mode  building  up  in  the  phonon  spectrum  i.e.  certain 
shear  waves  develop  a  very  weak  spring  constant  associated  with  the 
forces  returning  them  to  their  equilibrium  position  after  displace¬ 
ment.  These  correspond  to  the  very  directions  the  atoms  will  move 
when  they  undergo  the  crystal  structure  change.  At  the  transition 
temperature  the  atomic  positions  are  at  the  same  energy  so  that 
even  a  relatively  weak  compressional  or  shear  wave  can  unbalance 
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this  energy  difference  and  cause  large  displacements  of  the  atoms. 

When  the  displacements  are  large,  anharmonicity  becomes  effective 
and  this  in  effect  enables  the  long  wave  length  acoustical  waves 
to  be  quickly  damped  out  and  converted  into  heat.  From  the  point 
of  view  of  damping  materials,  the  conversion  of  vibrations  into 
heat  is  very  desirable  since  the  energy  of  a  phonon  is  proportional 
to  its  frequency  so  that  one  thermal  phonon  (1013  hertz)  can  absorb 
1010  acoustical  phonons  '103  hertz).  The  conversion  of  relatively 
large  quantities  of  vibrational  energy  produces  very  small  temperature 
rises. 


Any  theoretical  treatment  of  anharmonicity  and  "soft 
modes"  in  real  crystals  is  prodigious  so  we  consider  it  useful  to 
examine  a  very  simple  system  just  to  see  how  anharmonicity  couples 
normal  modes  together  and  affects  certain  physical  properties. 

The  system  we  have  chosen  is  that  of  the  one-dimensional 
coupled  oscillators  each  of  identical  mass  M  coupled  between  two 
rigid  (or  non-rigid)  walls  by  three  identical  springs.  The  potential 
energy  of  each  spring  is  represented  by  a  quadratic  term  (harmonic) 
and  two  anharmonic  terms,  one  cubic  (stronger  in  compression)  and 
one  quartic  so  that  the  total  potential  energy  of  the  coupled  oscilla¬ 
tors  becomes  (rigid  walls) 

V  *  KXj2-b3^2x^3  +g2Xi'<  +  KX22  b3^2X23  +g2x24  +  ■£•  k(xi-X2)2 

+b3^2(xi  -  x2)3  +  g2(xj  -  X2)1*  (1) 

where  ic ,  b  and  g  represent  the  spring  constants  respectively  of  the 
quadratic,  cubic  and  quartic  terms.  Since  we  wish  to  examine  the 
effect  of  the  anharmonic  terms  on  the  normal  modes  which  exist  in 
the  absence  of  the  cubic  and  quartic  terms  of  equation  (1)  we  trans¬ 
form  to  the  system  of  normal  coordiantes  ?  and  n  given  by 

5  -  +  x2)  ^ 

n  *  i(xi  -  x2)  J  (2) 

so  that  the  potential  is  now 

v  *  k(52  +  3n2)  -  6b3^2n(C2  -  n2)  +  2g2(£2  +  3n2)2.  (3) 

quadratic  cubic  quartic 

The  second  and  third  terms  on  the  right-hand  side  of  equation  (3) 
representing  the  cubic  and  quartic  terms  will  be  considered  small 
relative  to  the  quadratic  tetui  so  that  we  can  use  a  perturbation 
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treatment . 


In  the  absence  of  the  perturbation,  the  Schrodinger 
equation  is  soluble  exactly,  leading  to  a  wave  function  which  is 
a  simple  product  of  wave  functions,  one  for  each  mode, 


(n) 

m  n 

h  (a£j)  exp(  -  az£2) 

<l> (n)  ■  (  i/3*-  ^  ^  H  (31/4an)  exp(  -  /3a2n2) 

n  \2nn!  /it  /  n 

a  =  /2ir(Mx)  *^2 


C4) 


where  H  (y)  is  the  Hermite  polynomial.  The  eigenvalues  of  the  system 
are 


E  *  (m  +  -jjr)hv  +  (n  +  -J-)/3hv 


harmonic  approximation 


where  the  quantum  numbers  m  and  n  assume  integral  values  including 
zero,  and  refer  respectively  to  the  £  and  n  normal  modes.  The  £ 
normal  mode  represents  simple  harmonic  motion  in  a  system  in  which 
the  observer  is  always  positioned  at  one-half  the  algebraic  sum 
of  the  the  real  displacements  x^  and  X2,  while  the  n  normal  mode 
Is  the  same  for  the  average  of  the  difference  of  the  displacements. 
The  natural  frequency  of  the  n  normal  mode  is  v3  times  the  natural 
frequency  v  of  the  5  normal  mode,  equation  (5). 


When  we  apply  the  perturbations  of  equation  (3)  we  find 
the  cubic  term  vanishes  in  first  order.  We  have,  therefore,  calcu¬ 
lated  the  cubic  term  to  second  order  and  the  quartic  term  to  first 
order.  The  new  set  of  eigenvalues  becomes 

E=(m  +  +)hv+(n  +  4-)/3hv-  ^Trr-CT-Am-Sm^O '7mn+3*2+15n2+25'4n) 

a  hv 

quadratic  term  cubic  term 

2 

+  fir(6-5m  +  3m2  +  6*9nn  +  7’7  +  9n2  +  12*5n).  (6) 

Si 

quartic  term 
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New  terme  now  appear  in  the  eigenvalues  which  are  products  of  the 
quantum  numbers  m  and  n  and  which  are  associated  respectively  with 
the  5  and  n  'normal  modes’ .  Thus  one  rerult  of  anharmonicity  is  that 
the  excitation  of  one  of  the  ’normal  modes’  affects  the  energy  levels 
of  the  other.  (The  words  ’normal  modes'  are  now  in  quotation  marks 
since  they  strictly  refer  to  the  harmonic  solutions) . 

We  can  select  values  of  b  and  g  to  be  approximately  those 
obtained  from  an  expansion  of  the  Morse  potential  since  this 
potential  provides  a  good  fit  to  diatomic  molecules 


V  -  D  fexp  {-2X  -2  exp  (-*  (|5)‘/2}] 


(7) 


Mcrse  potential 


where  D  is  the  depth  of  the  potential  well  and  k  the  spring  constant. 
Expanding  around  the  minimum  one  obtains 


T  [«2  - 


.3/2x3 
T/T  + 


14D  J 


(2D) 1 ' 2  1  24D 

which  by  a  term-by-term  comparison  with  equation  (1)  gives 


b3/2 


g" 


k2/2 

WP77 

7kz 

48D 


(8) 


so  that  for  hv/D=10  2  t  equation  (6)  becomes 

E  =  hv[(m  f  £)+/3(n  +  £)+  10~2(0-02m  +  0-81m2  -  l*58mn 

+  0-69  -  0*57n2  -  0-57n)] 

energy  levels  for  a  Morse  potential. 

t  In  a  typical  case  hv  is  of  order  k9  (where  9  is  the  Debye  tempera¬ 
ture)  and  D  is  of  the  order  of  the  cohesive  energy,  so  that  a  valve 
of  10  2  for  hv/D  's  reasonable. 
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While  the  levels  of  low  quantum  numbers  are  only  slightly  altered,  the 
higher  levels  are  altered  significantly  owing  to  anharraonicity,  and 
in  some  cases  levels  even  cross.  It  is  clear  from  equation  (9)  that 
the  energy  levels  in  the  anharmonic  treatment  can  not  be  simply  re¬ 
presented  by  a  modified  harmonic  spectrum,  i.e.  one  in  which  the  two 
natural  frequencies  v  and  /3v  are  modified  with  increasing  tempera¬ 
ture  so  as  to  preserve  the  simple  form  of  equation  (5) .  This  must 
be  considered  an  artifice  that  avoids  the  basic  aspect  of  anharmon- 
icity  (i.e.  a  change  in  the  spectrum  of  eigenvalues). 

We  may  now  employ  thermodynamics  to  evaluate  the  high 
temperature  specific  heat  for  the  coupled  oscillators,  equation 
(6).  Setting 

b3/a6(hv)2  *  c  ^  1 

g2/a4(hv)  *  q  ^  1 

and  recalling  that 

.  -x 

9  -mx  1  +  e 

re  *  e  - 

(1  -  e“x)3 

we  obtain  the  partition  function  f  at  high  temperatures,  i.e. 
kT  }  hv  (neglecting  zero  point  terms  which  do  not  contribute  to 
the  specific  heat): 

£  -  7§£T  [l  +  22c  -  13  -7q  +  ^.  - 
X  *  hv/kT 

so  that  the  specific  heat  becomes  at  constant  length 

CL  “  kT  f^“(T  ln  f)*2k  (X  +  (l4c  "  10q)  hT-}  * 

Thus  in  agreement  with  the  classical  treatment  the  leading  term  in 
the  departure  of  the  specific  heat  from  the  Dulong-Petit  value 
of  2k  is  linear  in  T  either  positive  or  negative  depending  on  the 
magnitude  of  the  cubic  and  quartic  terms  in  the  potential.  For 
a  Morse  type  potential  (c/q  *6/7)  the  term  is  positive. 

If  we  allow  the  walls  to  move  so  that  the  length  of  each 
spring  is  charged  by  6 ({positive  for  compression),  new  terms  appear 


(ID 


(12) 


(cubic) 

(quartic) 


(10) 
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in  the  potential  which  are  added  to  those  in  equation  (3) .  The 
leading  new  terms  are 

k62  +  (66b3/2  +  12g262)(52  +  3q2)  (13) 

so  that  by  replacing  k->-k  +  66b3/2  +  12g262  for  the  £  and  q  modes  and 
adding  -|k62  to  equation  (6)  we  have  the  new  set  of  eigenvalues  as  a 
function  of  6 .  If  we  evaluate  the  high  temperature  free  energy  for 
this  set  of  eigenvalues  and  minimize  the  free  energy  relative  to  6, 
we  obtain  in  the  high  temperature  region 

2b3/2kT 

6  i?  (14) 


showing  that  the  system  expands  with  temperature.  The  two  frequencies 
of  the  system  v  an;<  /;,  are  now  temperature  dependent  and  become 


/,  6b3/2kT  \ 

v  * v  V 1  ■  -?—) 

/3v  »  /3v  (l  -  — ^-2k-Tj 


£  mode  \ 

f* 

n  mode  I 


For  a  soft  mode,  to  be  effective  the  cubic  and  quartic  terms 
must  become  proportionately  greater  than  the  quadratic  term. 
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